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1 Introduction 



An important lesson of nonlinear dynamical systems theory is that solutions originating 
from different regions of state space can (and often do) have qualitatively different modes 
of behaviour. In this sense, the initial conditions become an important ingredient in de- 
termining the dynamical evolution of the system. This then raises the question of how 
one determines the set of initial conditions that gives rise to a particular dynamical mode 
of behaviour. The answer to this question requires a full understanding of the underlying 
dynamics which, in the case of general relativity or string theory, is extremely difficult to 
establish, given the high degree of nonlinearity and complexity involved in these theories. 
There is also the added difficulty that despite recent progress in our understanding of M- 
theory, a definitive non-perturbative formulation of quantum gravity still remains to be 
developed (see, e.g., Ref. [0 for recent reviews). 

This question is particularly relevant in cosmological studies of string or M-theory. 
The low-energy effective action of the Neveu-Schwarz/Neveu-Schwarz (NS-NS) bosonic 
sector of string theory contains a multiplet of massless fields {g^u, <t>, B^ v } ||. The vacuum 
expectation value of the dilaton, (ft, determines the string (gravitational) coupling, g 2 s = , 
the gravitational field is determined by the metric, g^, and is an antisymmetric, two- 
form potential. To date, attention has focused on the spatially homogeneous, orthogonal 
models, where the dilaton field is constant on the surfaces of homogeneity ||, || . However, 
these models apply on scales just below the string scale, and it is precisely this region where 
spatial inhomogeneities may be important. 

A study of inhomogeneous string cosmologies is therefore necessary if further progress 
is to be made in addressing the question of whether our universe arose out of generic initial 
conditions. The purpose of the present paper is to derive and study a wide class of 'Einstein- 
Rosen' string cosmologies with non-trivial dilaton and two-form potential. Einstein-Rosen 
metrics are interesting for a number of reasons ||. (For reviews see, e.g., @,|). Spatial 
homogeneity is broken along one direction and they admit an Abelian isometry group, G2, 
that acts on two-dimensional spacelike orbits. They represent a natural generalization 
of the spatially homogeneous Bianchi cosmologies |], 1Q]. Density perturbations in the 
early universe can be analysed with these backgrounds and the propagation and collision 
of gravitational plane waves on homogeneous space-times can also be studied within this 
context [11, [l^,[l3|, 14]. Finally, it has been conjectured that G2 metrics represent a leading- 
order approximation to more general solutions near the singularity [15]. 

There has also been an interesting recent development within string cosmology, namely 
the pre~big bang scenario |16| , according to which the rapid increase of the string coupling 
drives an accelerated, inflationary expansion. The central postulate of this scenario is that 
the initial state of the universe is in the perturbative regime of small coupling and curvature. 
This leads to an inflationary phase for sufficiently homogeneous initial conditions [17, 15]. At 
present, the question of whether in general large spatial inhomogeneities have a significant 
effect on the naturalness of such initial data is unresolved. 

Recently, Barrow and Kunze [19] studied a class of inhomogeneous generalizations of 
Bianchi I string cosmologies and Feinstein, Lazkoz and Vazquez-Mozo |2(| derived an inho- 
mogeneous model by applying duality transformations on the locally rotationally symmetric 
(LRS) Bianchi type IX cosmology. In this paper we consider the G2 inhomogeneous general- 
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izations of the Bianchi type VI/,, universe, where h < is the group parameter. This Bianchi 
model is interesting because it has a non-zero measure in the space of homogeneous initial 
data and includes the Bianchi type III as a special case (h = —1) |2l]]. Furthermore, the 
most general spatially homogeneous solutions of the (one-loop) string equations of motion 
are the Bianchi types III and VI/,, where h = {0, —1/2, —2} [[H]]. It can be shown that these 
models contain the maximum number of eight free parameters. 

We employ non-compact, global symmetries of the field equations to generate inho- 
mogeneous solutions with a non-trivial two-form potential. When the metric admits two 
commuting spacelike Killing vectors, there exists an infinite-dimensional symmetry on the 
space of solutions that may be identified infinitesimally with the 0(2,2) current algebra 
[22, 23]. This symmetry reduces to the Geroch group, corresponding to the SL(2,i?) cur- 
rent algebra, when the dilaton and two-form potential are trivial The global SL(2, R) 
'S-duality' [25] and 0(2, 2; R) 'T-duality' [|(| are contained within this symmetry[|. Appli- 
cation of both these symmetries leads to new, inequivalent solutions. 

This paper is organised as follows. In Section 2, we derive inhomogeneous G2 string 
cosmologies from a general class of Einstein-Rosen models where the two-form potential is 
trivial. Two families of solutions representing inhomogeneous generalizations of the Bianchi 
type Vlft universe are found in Section 3 by directly solving the field equations. The 
asymptotic behaviour of such models is studied in Section 4. We conclude with a discussion 
in Section 5. 



2 Einstein— Rosen String Cosmology 



2.1 String Effective Action 

Fundamental strings sweep geodesic surfaces with respect to the string-frame metric, 
The four-dimensional, string effective action for the NS-NS fields is given by 



S 



d X\ 



-ge 



R + (V<A) 2 - ^H^ X H^ X 



(2.1) 



where H^ u \ = d^B uX ] is the field strength of the two-form potential. In general, the 
effective action will also include moduli and vector fields arising from the compactification 
from higher dimensions. The action is also expected to include a potential term, V(4>), 
arising from the non-perturbative sector of the theory, however, the form of such a potential 
is as yet unknown. These additional terms are neglected in what follows. 

In order to take advantage of the highly developed framework of general relativity and 
its many known exact solutions, it is often more convenient to work in the Einstein frame, 
where the dilaton field is minimally coupled to gravity. This is achieved by making the 
conformal transformation 

g^u = eT^g^u- (2.2) 

Action ([O]) then takes the form 



S 



R 



(2.3) 



We refer to the groups SL(2, R) and 0(2,2;_R) as the S- and T-duality groups, respectively, although 
at the non-perturbative level the dualities are the discrete subgroups SL(2, Z) and 0(2, 2; Z). 
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In four dimensions the field strength of the two-form potential is dual to a one-form: 

H» vX = e^ AK e 2 ^V K <7, (2.4) 

where e^ uXK is the covariantly constant four-form J25[ |. The field equations can then be 
derived from the dual action 

(2.5) 



where a may be interpreted as a pseudo-scalar 'axion' field. 

The generalized Einstein-Rosen G2 metric is defined in the Einstein frame by the line 
element g [7| 

ds 2 e = h a/3 (x' y )dx a dx p + % b (x' y )dx a dx b , (2.6) 

where all components are taken to be independent of the spatial coordinates x a = (x,y). 
The two commuting, spacelike Killing vectors are d/dx and d/dy and h a p represents the 
longitudinal component of the metric. The metric on the surfaces of transitivity is denoted 
by Jab and the gradient = (det^) 1 ^ 2 determines the local behaviour of the model. 
Solutions represent cylindrical and plane gravitational waves if is globally spacelike or 
null, respectively [0, 11 1. Cosmological models arise when is timelike or when the sign 
of K^K^ changes || . 

The longitudinal metric is conformally flat and the line element may therefore be written 
in the form 

ds 2 e = e s (-df + dz 2 ) + % b dx a dx b , (2.7) 

where / = / (£, z) determines the longitudinal part of the gravitational field. The corre- 
sponding line element in the string frame is given by 

ds 2 s = e^ +f i-di 2 + dz 2 ) + T ab dx a dx b , (2.8) 



where the transverse metric T ab = e^j ab has determinant T = detr a fe. We assume through- 
out this work that all massless fields are independent of the coordinates x a . Thus, the 
metric (2^) also represents a G2 model. 

A considerable simplification occurs in the field equations when the transverse metric is 
diagonal and separable. In this case, the Einstein and string frame metrics may be written 
in the form 

ds 2 e = e f {-df + dz 2 ) + £{e p dx 2 + e~ p dy 2 ) (2.9) 

and 

ds 2 s = e f+ ^ (-^ 2 + dz 2 ) + ^ (e p dx 2 + e~ v dy 2 ) , (2.10) 

respectively, where p = p(£, z) represents the transverse part of the gravitational field. In 
this case, the volume of the transverse space in the string frame is determined by 

T = e 2 e 2 *. (2.11) 

In some settings, it proves convenient to define new variables 

z = e~ 2Z cosh(2t) 

£ = e~ 2Z sinh(2t), (2.12) 
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which transform the metric ( |2,9| ) to 

ds 2 e = Ae h {-dt 2 + dZ 2 ) + e~ 2Z smh(2t)(e p dx 2 + e~ p dy 2 ), (2.13) 
where h = f — 4Z. 

In the present context, an inhomogeneous string cosmology is parametrized by the mass- 
less degrees of freedom {g^ u ,(p,B^}. A vacuum solution to Einstein gravity may then be 
represented by {g^, 0, 0} and a solution with a trivial two-form potential by {g^ u , 4>, 0}. We 
refer to this latter class of solution as 'dilaton-vacuum' solutions. In general, the one-loop 



field equations of motion for G2 backgrounds derived from the actions fl2,l| ) and (2.5) are 
difficult to solve. In view of this, we employ the non-compact, global symmetries that arise 
when the metric admits two Abelian isometries to generate a wide class of inhomogeneous 
string cosmologies with a non-trivial two-form potential from dilaton-vacuum solutions. 

2.2 0(2,2) Symmetry 

The global 0(2, 2) symmetry applies when there exist two Abelian isometries and the only 
non-trivial component of the two-form potential is B xy = B xy (^, z) pq| . This symmetry is 
manifest in the string frame and generates fractional linear transformations on the two-form 
potential and the components of the transverse metric, T a i,. When the transverse metric is 
non-diagonal, the four degrees of freedom {B xy , T a b} parametrize the 0(2, 2)/[0(2) x 0(2)] 
coset [26]. The isomorphism 0(2,2) = SL(2,i2) x SL(2,i?) then implies that these may be 



arranged in terms of two complex coordinates |2 

+ I — 

r r 

1 yy L yy 



t^ + z^ (2.14) 
1 yy L yy 

p = B xy + iVf. (2.15) 



Suppose the metric ( 2.1C| ) represents a vacuum or dilaton-vacuum solution for some 



appropriate form of the dilaton, (j) = 4>(£,i z )- An 0(2,2) transformation is then generated 
in terms of the two SL(2, R) transformations: 

r = r (2.16) 
p = p, (2.17) 





ap + b 


p = 


cp + d' 




o!t + b' 


f = 


c'r + d r 


b'c' = 1. 


Under a 



transforms to ^ 

S> = (2.18) 

and the longitudinal part of the string frame metric remains invariant, i.e., 

/=/ + 0-£ (2.19) 
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Since the transformations ( p.!7|) leave the two-form potential invariant, (2.16) must be 
employed to generate such a field from a dilaton-vacuum seed solution. Applying ( 2.16 ) 
implies that 

f = (2.20) 

(d 2 + c 2 r) 2 

- acT + bd . , 

Bxy ~ ^FT^' ( 2 - 21 ) 

(2.22) 



d 2 + c 2 r , 



where the last transformation follows from fl2,ll| ) and ( 2.20|) . The dual metrics in the string 
and Einstein frames are then given by 

ds* = e?+* (-df + dz 2 ) + d2 (e?dx 2 + e~Uy 2 ) (2.23) 



and 



dsl = e f (d 2 + c 2 e^) (-dC 2 + dz 2 ) + £ (e^dx 2 + e^dy 2 ) , (2.24) 

respectively. We remark that the transverse metric in the Einstein frame, 7^, is invariant 
under the transformations ( 2.16| ). When d = and c = 1, the volume of the transverse 



space in the string frame, as given by Eq. ( 2.2C| ), is inverted. The transformations ( 2.16 ) 
therefore represent a 'T-duality'. 

2.3 SL(2,R) Symmetry 

The global SL(2, R) symmetry of the string effective action ( |2.1| ) becomes manifest in the 
Einstein frame. The action (|2.5|) may be written as a non-linear sigma-model, where the 
dilaton and axion fields parametrize the SL(2, R)/\J(1) coset [25]. The effective action is 



therefore invariant under global SL(2, R) transformations. These act non-linearly on the 
complex scalar field \ = a + such that the transformed field is given by 

_ A X + B 

X = C^D' (2 - 25) 

where {A, B,C, D} are real numbers satisfying AD — BC = 1. The Einstein frame metric 
transforms as a singlet under this SL(2, R) transformation and the dual string frame metric 
is therefore given by 

ds 2 = e^ds 2 . (2.26) 

For the special case where C 2 = 1 and a = —D/C, the SL(2,i?) transformation fl2,25| ) 
yields 4> = which corresponds to an inversion of the string coupling g s = gj 1 . The 
transformations ( |2.25| ) therefore represent a strong/weak-coupling 'S-duality'. 

Starting with a = 0, a solution with non-trivial axion field may be generated directly 
from a given dilaton-vacuum solution of the generic form ( |2.9D by applying Eq. ( p. 25 ). The 
dual solutions are given by 

ds 2 = e^+f (-de 2 + dz 2 ) + (fdx 2 + e~ p dy 2 ) (2.27) 
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= C 2 e _ * + I>V (2.28) 
ACe~* + BDe* 
C 2 e-^ + D 2 e<t> 

To summarize thus far, we have seen how two inequivalent classes of inhomogeneous 
string cosmologies with a non-trivial two-form potential can be derived from a given 
dilaton-vacuum solution, by employing the non-compact, global SL(2, R) and 0(2,2) sym- 
metries of the model. The dual solutions may be referred to as 'dilaton-axion' cosmologies. 
In all cases, they are parametrized in terms of the functions {f,p, 4>} that define the seed 
dilaton-vacuum solutions. 

Thus, the asymptotic behaviour of these models can be investigated directly once the 
seed solution has been specified. We therefore proceed in the following Section to derive two 
classes of inhomogeneous dilaton-vacuum cosmologies that may be viewed as generalizations 
of the homogeneous Bianchi type VI/i universe. 

3 Inhomogeneous Dilaton— Vacuum Cosmology 
3.1 Cosmological Field Equations 

When the two-form potential vanishes, action (^1]) reduces to that for a massless, minimally 
coupled scalar field. For the metric Q2.9|) the field equations then take the form pq] 



f = -i + Up 2 +P ,2 + 4> 2 + ^) (3.1) 



2£ 2 

/' = t (pp' + #') , (3.2) 

P+]p-p" = 0, (3.3) 

+ ~J> - 4>" = 0, (3.4) 

in which overdots (primes) denote differentiation with respect to the timelike variable £ 
(spacelike variable z). Eqs. fl3.3| ) an d (|3.4[ ) are the integrability conditions for the system 
(|3.l| ) and ([^). The field equations (]3.l|) -( |3.4|) are invariant under the simultaneous inter- 
change p <-> <f). Indeed, the wave equations ( |3.3| ) and Q3.4| ) are formally equivalent to the 
cylindrically symmetric wave equation in flat space, for which the general solution is known. 
A number of cosmologically relevant solutions satisfying different boundary conditions have 
been considered previously ^, ^9[. An important feature of these equations is that they 
are linear, which implies that new solutions may be constructed from superpositions of 
known solutions. 

Here we consider a linear superposition of solutions of the form: 

p = Hn(£) - m cosh" 1 (~J + e x J™ [ Cl (z)J (lO + c 2 (z)N (lO] dl (3.5) 
</) = aln(0-/?cosh- 1 (£) +e 2 J™{c 3 (z)J (lti) + c 4 (z)No(lO] dl, (3.6) 
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where {a, /3,k,m} and {£1,62} are constants, Jo and Nq are zero-order Bessel functions of 
the first and second kind and the coefficients q = Cj(z) are defined by 



d(z) = Qcos{lz) + Disin(lz) (3.7) 

c 2 (z) = F t cos(Zz) + Gi sin(/z) (3.8) 

cs(z) = Hi cos(lz) + L[ sin(Zz) (3-9) 

c 4 ( z ) = Uicos(lz) + Vism(lz), (3.10) 

where Ci, Di, etc., are arbitrary constants for each value of I. In both Eqs. ( |3.5[) and ( |3.6D , 
the last term represents the most general separable solution to Eqs. ( |3.3| ) and (p.4|). The 
remaining two terms represent other, in general inhomogeneous, solutions^. The importance 
of these other solutions is that they include a number of spatially homogeneous Bianchi 
models as special cases. The Bianchi models admit a three-dimensional group of isometries, 
G3, that acts simply transitively on three-dimensional spacelike orbits The G3 contains 
an Abelian subgroup G2 for the types I-VII/i and the LRS types VIII and IX || ^Cj] . Certain 
G2 models may therefore be viewed as inhomogeneous generalizations of these Bianchi 
cosmologies. In particular, for solutions of the form Q3,5| )-( |3~l~0| ), one has the following 
sub-classes of solutions: 



1. The class with ei = e2 = m = (3 = corresponds to the homogeneous orthogonal 



Bianchi type I models containing a stiff perfect fluid [31]. They reduce to the Kasner 
solution in the vacuum limit (a = 0) p^j. 

2. The class with m = (3 = 0, e\ = 62 = 1 corresponds to an inhomogeneous generali- 
sation of the Bianchi I models. Charach and Malin first considered solutions similar 



to these by imposing a three-torus topology on the spatial sections [28]. This effec 



tively converts the integral over I into an infinite sum, which also has the consequence 



of simplifying the solution of the remaining field equations for /. Adams et al. [12] 
further considered the vacuum case (0 = 0). 

3. The class with ei = €2 = 0, subject to the constraint 1 — a 2 + m 2 — k 2 — 3 = 0, where 
k = (— h)~ l l 2 . In general, these models represent tilted stiff perfect fluid Bianchi type 



VI/i cosmologies [33], since the fluid velocity vector associated with the dilaton field 
is not orthogonal to the group orbits (surfaces of homogeneity). They reduce to the 
Bianchi type III and V models when k 2 = 1 and k = 0, respectively In the 

vacuum limit (a = /3 = 0), the solution reduces to the Ellis-MacCallum type VI^ 
cosmology p5[. 



4. The class with e\ = €2 = corresponds to the inhomogeneous generalisations of the 
Bianchi type III, V and VT/j models first considered by Wainwright et al. [33]. In 



general, however, these models suffer from spacelike singularities and their status as 
physical cosmologies is uncertain M. 



2 The first term is also present in the third term and is therefore not a different solution. It is separated 
out for later consideration. 
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An alternative way of considering inhomogeneous generalisations of the Bianchi type 
III, V and Vlft models is to specify e± = €2 = 1 and impose the constraint 



(3 2 - a 2 + m 2 - k 2 - 3 = 0. (3.11) 



Vacuum solutions of this type were considered by Adams et al. [13], who concluded that 
the inhomogeneous structure of the initial cosmic singularity could evolve into gravitational 
waves propagating over a homogeneous background at late times. In view of the ambiguities 



associated with interpreting the Wainwright et al. [33] solutions in a cosmological context, 
our primary interest in the present paper is in this new class of models. We also remark 
that a subset of inhomogeneous Bianchi I models (item 2) is also included within this class 
since these latter models correspond to the particular solution (3 = m = 0. 

The remaining field equations ( |3.1| ) and ( |3.2| ) for the longitudinal metric function / 



may now be solved in principle by substituting in the derivatives of Eqs. (BJj) and (pM 



However, solving Eq. (3.1) is non-trivial, because the right-hand-side contains integrals over 
I that originate from the integral wave-train terms of Eqs. ([3j]) and (|3.6| ). Unfortunately, 
these can not be expressed in a closed form and some simplifying assumptions must therefore 
be made in order to proceed analytically. 

In view of this, we consider two separate schemes. In the first, we restrict the analysis to 
a single mode, I, which could be viewed as dominating over all the other modes. The choice 
of one mode allows the integrals over I to be dropp ed^} In the second case, we assume the 
amplitudes of each of the modes in Eqs. (|3.5|) and (|3.6[) are equal, i.e., we specify C% = C, 
Di = D, etc. 



We first integrated Eq. (3T) to deduce an expression for / containing an unknown 
function of integration fi(z). In each of the cases we considered, it was found that Eq. 
was then trivially satisfied for constant fi(z). We now proceed to present the two 



classes of solutions, together with their homogeneous limits. 



3.2 Homogeneous Solutions 

The homogeneous limit of these solutions is determined by specifying ei = €2 = in Eqs. 



and (|3.6|) |33| . Eq. (|3.l| ) can then be integrated to yield the longitudinal component 
of the gravitational field in the metric (|2.9[): 



hom 



Cl + l(a 2 + (3 2 + k 2 + m 2 - 1 ln(£) --((3 2 + m 2 ) ln(z 2 - £ 2 ) 



(a/3 + km) cosh 



-1 



(3.12) 



where C\ is an arbitrary constant of integration and the constraint equation ( 3.1 1| ) applies. 
In terms of the variables ( [2.12 ) this component is given by 



hho 



d + - (a 2 + (3 2 + k 2 + m 2 - l) In sinh(2t) + (km + a/3) In tanh t 



(3.13) 



3 Particular solutions of the field equations that consist of superpositions of two o r m ore m odes may also 
be found by employing a discrete summation over I in place of the integral in Eqs. (3.5) and (3.6). 
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in Eq. ( p. 13 ), The transverse component of the metric and the dilaton field respectively 
take the following forms 



Phom = —2kZ + k\nsmh(2t) +mlntanht (3-14) 

<fthom = — 2ctZ + alnsinh(2t) + /^lntanhi. (3.15) 

The hypersurfaces t = constant represent the surfaces of homogeneity. Since for a / 0, 
the dilaton field depends on the spatial variable, Z, the fluid velocity as measured by 
= (j> u/(— 4',v c t >, ' J ) 1 ^ 2 is n °t orthogonal to the group orbits. Thus, this solution may be 
interpreted clS 8l tilted stiff fluid type VI^ solution. 



3.3 Single-mode solutions 

After restricting the analysis to a single mode, I, the integral of Eq. ( |3.1[ ) is readily found. 
The resulting solution may be expressed in terms of its homogeneous, gravitational wave 
and scalar wave components: 

/ = fhom ~l" fgw ~\~ fswi (3.16) 

where the homogeneous component, fhom, is given by Eq. ( 3.12| ), the gravitational wave 
component is given by 

-2 



fgw = +j(KJo(/0 + c 2 iVo(^)] 2 + [c / 1 J 1 (/0 + c / 2 iVi(/0] : 



I 2 

+ 



e 2 (Mo(Z£) + c 2 N (lC)] 2 + [eMIO + c 2 N 1 (lO\ 2 



4 

~Z [ciMlO + c 2 N (lO] [cMlO + 02^(10} + k ( Cl J (Z£) + c 2 N (lO) 

J Vz 2 - e J y/^e 



-mlzci 



and the scalar-wave component f sw is determined by making the substitutions 

k — > a, m — > 0, c\(z) — ► 03(2), c 2 {z) — » c±(z) (3.18) 

in Eq. (3.17). The gravitational and scalar wave components encode all the inhomogeneous 
contributions to the longitudinal component of the gravitational field. The integrals in 
the expressions for f gw and f sw can actually be performed, provided the integrand is first 
expressed as a series. The result is a series which reduces to a closed form expression in the 
asymptotic late-time limit. 



3.4 Equal-amplitude solutions 

When all the modes in the integral wave-trains of Eqs. ( |3,5j ) and ( |3.6| ) have equal weighting, 
we may employ the identities 



1 



if f(lz) = sin(Zz) 



MiOf(iz)di = { V^e (3.19) 

[ 0, if f(lz) = cos(Zz) 
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N (lOf(lz)dl 



2 cosh 1 (z/£) 



, if f(lz) = sin(lz) 
if f(lz) = cos(lz) 



(3.20) 



in order to integrate Eq. ( |3.lD . To avoid unnecessary complications, we further restrict 
our attention to the class of solutions where the contributions from the Neumann functions 
vanish (c2 = 04 = 0)0. The transverse component of the gravitational field and the dilaton 
are then given by 



p = k m(£) — m cosh 1 ( - ) + 



F 



aln(£) - (5 cosh' 1 ( - ) + 



M 



(3.21) 



(3.22) 



respectively, where F and M are arbitrary constants. The longitudinal component of the 
gravitational field, /, may now be determined by substituting Eqs. ( 3.21 ) and ( |3.22| ) into 
Eqs. ([O]) and ( |3.2[) and integrating. One finds that 



_ , F 2 + M 2 t2 mF + pM kF + aM 

I - ham + 4(z2 _ ^ + z2 _ e Z+ ■ 13.23) 

In terms of the variables defined in Eq. ( 2.12Q , it is given by 

h = h hom + (kF + aM)e 2Z + (mF + (3M)e 2Z cosh(2t) + ^(F 2 + M 2 )e AZ sinh 2 (2t), (3.24) 

where hh om is given by Eq. ( 3.13j ). The transverse component and dilaton are given by 

p = -2kZ + kin sinh(2t) + m In tanh t + Fe 2Z (3.25) 

4> = -2aZ + a In sinh(2t) + f3 In tanh t + Me 2Z (3.26) 
when expressed in terms of these variables. 



4 Asymptotic Behaviour 



In this section we consider some aspects of the asymptotic behaviour of the inhomogeneous 



string cosmologies derived above. The axion field in the dual solution ( 2.27 )— ( 2.29 ) tends 
to a constant value in the limits <j) — > ±00. It is therefore dynamically negligible in these 
limits. Eq. ( 2.28| ) implies that for all solutions generated by the SL(2,i2) transformation 
( 2.25| ), there exists a lower bound on the value of the dilaton field if C and D are non-zero. 
This implies the existence of a lower (non-vanishing) bound on the string coupling which, 
in the context of M-theory, in turn implies the existence of a lower bound on the radius of 
the eleventh dimension^, IZn, since TZn oc e 



0/3 



d^2iD 



that the two-form 
00, i.e., when the volume 



In the class of dual solutions ( |2.23| ), it follows from Eq. 
potential tends to a constant in the limits where r — > and T - 
of the transverse space measured in the string frame becomes vanishingly small or arbitrarily 

4 The solution may also be found for the case when the Neumann functions are included. 
5 We are assuming implicitly that the extra six spatial dimensions are fixed. 
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large. Thus, the two-form potential effectively decouples from the field equations in these 
limits. The limiting behaviour of the dilaton field in the dual solution follows from Eq. 
( 2.22j ). In particular, when r — ► and d 7^ 0, the dilaton asymptotically tends to the form 
it has in the seed solution (up to a constant linear shift). The same conclusion applies for T. 
In this limit, therefore, the dual cosmology asymptotes to the original seed, dilaton-vacuum 
solution. In the opposite limit, where T diverges, the dual cosmology asymptotes to the 
solution that is derived by specifying d = in Eq. ( 2.16|) . 

The above discussion applies to any class of string cosmologies derived with the sym- 
metry transformations discussed in Section 2. To illustrate this further, we consider the 
asymptotic behaviour of the string cosmologies ( |2.23| ) derived from the equal-amplitude 
dilaton-vacuum, seed solution ( 3.21 )-( 3.23 ). 

We recall that in determining the asymptotic behaviour of cosmological models the 
choice of time gauge is important. Since the models considered here may be viewed as de- 
scribing inhomogeneous waves propagating over homogeneous Bianchi backgrounds, a rea- 
sonable measure of early and late times is provided by t in the coordinate chart {t, Z,x,y}. 
As we shall see, the asymptotic behaviour in the limits of small and large t corresponds to 
either r — » or 00 in the solution ( 2.23j ). This implies that the class of string cosmologies 
( 2.23| ) asymptotes between two dilaton-vacuum solutions, where the two-form is dynami- 
cally negligible. In effect, this field induces the transition between the two dilaton-vacuum 
limits. A similar conclusion holds for the dual solution ( [2. 27 ). 

In the limit, t — > +00, the relevant terms in the asymptotic forms of Eqs. ( |3,21| ), J3.22 ) 
and (3.24) are 



h 



p ps 2k(t - Z) + F e 2Z 
<f> ps 2a(t - Z) + M e 2Z 

( a 2 + f3 2 + k 2 + m 2_ 1)t + 1 (F 2 + M 2 )e MZ+t)_ 

16 



In terms of the coordinate pair {£, z}, this implies that 

F 



P 



feln(0 + 



M 

Vz 2 - e 



1 



/ « + + k z + m z - 1) ln(0 + 



(F 2 + M 2 )i 2 



(4.1) 
(4.2) 

(4.3) 



(4.4) 
(4.5) 

(4.6) 



2' ' r _/ ^ v " ' Is:- ef ' 

We note that, in view of Eq. ( [2.12 ), the singularity at z = £ only corresponds to (£, z) — > 
(00, cxd) and therefore cannot be reached in a finite t-time along any curve. 

In the limit t — > 0, the corresponding limiting forms of the metric components and 
dilaton field are 



p PS -2k Z + (k + m)lnt + Fe' 
<p ps -2aZ + (a + /3)lnt + Me 



(4.7) 
(4.8) 



2 



a 2 + (5 2 + k 2 + m 2 + 2{km + a/3) 



ln(i) + [(k + m)F + (a + 0)M] e 2Z . (4.9) 



11 



The limits of the determinant T are deduced by substituting Eq. ( |2.12| ) and Eq. (4.2) or 
into Eq. (|2.11|) : 



T(t +00) 
F(t -+ 0) 



exp 



A(a + l)(t - Z) + 2Me 



2Z 



4exp -4(1 + a)Z + 2Me 2Z + 2(1 + a + (3)lat 



(4.10) 
(4.11) 



The behaviour of the dilaton field in the limit where the determinant diverges is most readily 
determined by substituting Eq. (2.11) into Eq. (|2.22j ). We find that (ft ~ — (ft — 2 ln£ — 21nc. 
This is interesting, because for finite Z, Eq. ( [4.2D implies that (ft oc aln£ for t — > 00. For 
q>1, therefore, (ft ~ —(ft and, in this sense, the strongly coupled limit of the dual solution 
may be viewed as the weakly coupled limit of the seed solution, and vice-versa. 

In taking the early or late time limits of an inhomogeneous cosmological model, one 
has to take account of both direction and time. Here, we have inhomogeneity in the in- 
direction only. The question that then arises is whether the Z-dependent terms in the above 
expressions dominate over the t-dependent ones. For any finite Z, the time-dependent terms 
eventually dominate as t — ► {0, 00}. A possible ambiguity in the limit arises, however, if we 
allow Z(t) — > 00 sufficiently fast relative to the i-dependent terms. 

It is helpful to consider the simple set of straight lines V = {Z = pt + k \p, k £ 1Z}, 
as probes with which to study the asymptotic behaviour^]. In this case, it follows from Eq. 
( |4T0| ) that 

r(t^n\~i exp[4(l + a)(t-Z)] if 
l(t-ooj~| ew[2Me 2 Z] if 

The sign of p determines the term that dominates in Eq. ( 4.12j ), with p > (p < 0) 
corresponding to motion in the positive (negative) Z-direction while p = corresponds to 
moving along Z = constant trajectories. Thus, for p < 0, T — > 00 for a > — 1, while r — > 00 
for M > and p > 0. On the other hand, the determinant becomes vanishingly small if 
a < -1 and p < or if M < and p > 0. 

In the early time limit, Z — > k = constant for the paths we are considering. Eq. (4.11) 
then implies that T oc t 2{l+a+ ^ . It follows that r^0ifa + /?>-landr^ooifa + /3< 
— 1. The limiting forms of the dilaton and the transverse and longitudinal components of 
the metric are determined by allowing Z — > k in Eqs. fl4.7D -( fQj|) . Transforming to the 
synchronous frame by defining 



P < 
p > 



(4.12) 



dt'e h ^' 2 



(4.13) 



then implies that, for both seed and dual solutions, the G2 line-element (2.6) in either the 
Einstein or string frames qualitatively takes the form 



ds 2 



-dr 2 + A^r^dx 2 + A 2 {n)T a2 dy 2 + A 3 (n)r a3 dz 2 



(4.14) 



where the constants cij can be expressed in terms of {k, m, a, (3} and Ai depend on k. This 
solution represents an inhomogeneous generalization of the Kasner-Belinskii-Khalatnikov 



(KBK) solution [38], in the sense that at each Z, parametrised here by the constant k, the 
universe describes a particular KBK solution. 

6 The precise form of these curves is not important, their utility derives from the fact that they may be 
employed to probe the three dynamically important cases Z(t) — > ±00, Z(t) = finite. 
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Thus, the determinant of the transverse space in the string frame vanishes or diverges 
in the early and late time limits. The precise behaviour depends on the constants that arise 
in the seed dilaton-vacuum solution. These constants are arbitrary modulo the constraint 
equation (3.11). In both limits, the two-form potential asymptotes to a constant value and 
the dilaton field tends to its original form when r —* 0. 



5 Conclusion and Discussion 

In this paper we have employed the global symmetries of the string effective action to 
derive two families of inhomogeneous string cosmologies from a generalized Einstein-Rosen 
metric admitting an Abelian group of isometries, (?2- The solutions were parametrized 
in terms of the metric functions of the seed solution. Thus, the qualitative behaviour of 
these models can be determined directly from the asymptotic form of the original metric. 
Inhomogeneous generalizations of the Bianchi type VI/j string cosmologies containing a 
non-trivial dilaton and two-form potential were derived. In general, these fields induce 
inhomogeneities that may be viewed as scalar and gravitational waves propagating over a 
homogeneous background. 

One of the main applications of Einstein-Rosen models arises because many of the non- 
perturbatively exact superstring backgrounds constructed from the gauged Wess-Zumino- 
Witten (WZW) models admit two abelian isometriesQ. Hence, given our current knowledge 
of string cosmology and conformal field theory, the G2 Einstein-Rosen cosmologies derived 
within the context of the low energy effective action represent a set of models that is closely 
related to many exact string solutions in four dimensions. In this sense, it is of primary 
importance to understand the dynamical behaviour of these solutions. In particular, one 
question that arises is whether the G2 solutions to the low energy string equations of motion 
discussed in this work asymptote towards known WZW models in the early- or late-time 
limits. 

Since our solutions are valid in the weak coupling regime, they should correctly represent 
the asymptotic states of perturbatively exact string cosmological models at future, or past, 
timelike infinity. On the other hand, in the strong coupling regime, i.e. near the big- 
bang singularity, one expects that higher-order corrections to the perturbative theory, as 
well as non-perturbative string effects, should become increasingly important. Hence, in 
this regime the qualitative behaviour of these solutions may deviate somewhat from that 
of solutions derived from the full M/string theory. However, there are reasons to believe 
that G2 solutions should nevertheless provide a generic description of cosmological models 
in the vicinity of a singularity. A major incentive for this comes from the long standing 
conjecture of Belinski and Khalatnikov | 38| , 39 , 4C ] . This states that on the approach to the 



cosmological singularity, the generalized Einstein-Rosen metrics may play the role of the 
leading-order approximation to the general solution of conventional Einstein gravity. It is 
therefore important to study the behaviour of these models in the t — > limit as well. 

In addition such solutions provide a useful framework within which to study a number 
of other topics in early universe string cosmology. In particular, they serve as a theoretical 



setting for investigating the pre-big bang inflationary scenario [16|. In the context of this 



7 See, for instance, Ref. |J7j f° r a review of the gauged WZW models. 
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scenario, one is usually interested in solutions where both the curvature and effective string 
coupling (the dilaton) diverge as t — > 0. So, from Eq. ( |4.8|) , it follows that one needs a+/5 < 
in order to satisfy the requirements for pre-big bang inflation in the string cosmologies 
generated from the equal-amplitude, dilaton-vacuum solutions. These solutions possess a 
Kasner-like early time limit which under time reversal makes them compatible with pre- 
big bang inflation. However, in order for the inflation to end, this scenario requires a 
mechanism for inducing a graceful exit into the standard, post-big bang phase. At present, 
this problem is unresolved and it is possible that the inhomogeneous singularity may pose 
a further complication in addressing the exit problem. It would be interesting to consider 
the homogenization of the universe within the context of these solutions. 

The initial conditions for the pre-big bang scenario are based on the assumption that 
both the curvature and coupling are sufficiently small enough to ensure that the universe 
is in the perturbative regime. An important question in connection with this scenario is 
the naturalness of these initial conditions [41]. An attempt has recently been made to 
address this question by conjecturing that the past attractor of the inflationary solution 
is likely to be that of the Milne universe [17]. If true, this would go some way towards 
establishing the 'naturalness' of the scenario. More general initial states were also recently 
considered [|^, |4^| . The solutions discussed above possess a range of possible late time 
limits, which in general are not Milne-like. This suggests that this class of models is not 
generally compatible with the conjecture of Ref. [17|. 

Furthermore, there exists a lower bound on the string coupling for models generated by 
the SL(2,J2) transformation ( p. 25 ) and this has implications for the range of initial values 
that such a parameter can take. This in turn leads to an upper limit on the amount of 
inflation that can occur before higher-order effects become significant [[4l|]. 

The solution generating techniques discussed in Section 2 can be incorporated into more 
general algorithms. Although the models presented in this paper break spatial homogeneity 
along one direction, they still exhibit a certain degree of symmetry and it is important to 
develop further techniques that lead to more general solutions. Recently, by extending a 
previous method [43], an algorithm was presented that generates inhomogeneous G\ scalar 
field cosmologies exhibiting a single isometry from matter filled and vacuum G2 backgrounds 
[44]. Such models break homogeneity in two spatial directions. The discussion of Ref. [44] 
was placed within the context of Einstein gravity and it would be interesting to adapt this 
algorithm to string cosmology. In principle, a family of inhomogeneous G\ string cosmologies 
could then be generated from the G2 solutions discussed in Section 2. 

The string effective action exhibits a further discrete symmetry when there exits a 
G2 isometry p2| . This 'mirror' symmetry interchanges the transverse metric degrees of 
freedom with the dilaton and axion fields and leads to a new solution with a different 



spacetime interpretation [22, [49]. When the axion field is trivial, the symmetry reduces 
to the simultaneous interchange p <-> 4> in Eqs. fl3,l[)-( |3~4|) . In particular, the negatively 
curved Friedman-Robertson- Walker string cosmology Q may be generated in this way from 
a vacuum Bianchi type V model [45]. A third class of inhomogeneous string solutions with a 
non-trivial two-form potential may therefore be found by applying this discrete symmetry 
to the Gi backgrounds derived in this work. 

The G2 backgrounds we have discussed are parametrized by non-trivial fields from the 
NS-NS sector of the string effective action. This sector is common to all five perturbative 
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string effective actions and the solutions may therefore be viewed as truncated solutions of 
both the heterotic and type II theories. The type IIB theory also has a non-trivial Ramond- 
Ramond (RR) sector, consisting of an additional axion field and a two-form potential ||. 
These fields differ from those of the NS-NS sector in that they do not couple directly to the 
dilaton field in the effective action |36|, . There are further symmetry transformations that 
can be applied to generate a non-trivial RR sector from a given NS-NS background p47[ |. Our 
solutions therefore represent seeds for investigating the role of RR fields in inhomogeneous 
string cosmologies. To date, such fields have only been studied in an homogeneous setting. 
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